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Abstract 

Interpopulation hybridization can lead to outbreeding depression within affected populations due 
to breakdown of coadapted gene complexes or heterosis in hybrid populations. One of the 
principal methods commonly used to estimate the level of developmental instability (DI) is 
fluctuating asymmetry (FA). We used three genetically differentiated Drosophila subobscura 
populations according to inversion polymorphism analysis and measured the variability of 
sternopleural bristle number and change in FA across generations P, Fl, and F2 between intra- 
and interpopulation hybrids of D. subobscura. The mean variability of sternopleural bristle 
number in intra- and interpopulation hybrids of D. subobscura across generations cannot 
determine whether the changes at the level of developmental homeostasis are due exclusively to 
genomic coadaptation or to heterozygosity. Phenotypic variance (V p ) and FA of sternopleural 
bristle number was higher in interpopulation than in intrapopulation hybrids across generations. 
Fl hybrids were more developmentally stable compared to each parental population in both intra- 
and interpopulation hybrids. The most probable mechanism providing developmental 
homeostasis is heterozygote or hybrid superiority, also called overdominace. However, V p was 
higher and FA lower in the F2 generation when compared to Fl, due mainly to crossing-over in 
the formation of F2. 
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Introduction 



Hybrid zones occur when two genetically 
distinct populations meet, mate, and produce 
offspring of mixed ancestry (Barton and 
Hewitt 1985; Ross and Robertson 1990). 
These zones can occur in disrupted habitats 
through interpopulation hybridization, where 
previously isolated populations become 
intermixed and hybridize. Anthropogenic 
activities in natural ecosystems may increase 
opportunities for interpopulation hybridization 
through habitat disruption. Hybridization 
between individuals from different 
populations can lead to heterosis, also called 
overdominace, in the first generation, 
followed by outbreeding depression in the 
consecutive generation due to the breakdown 
of coadapted gene complexes (Dobzhansky 
1950; Andersen et al. 2002; Edmands 2007). 
There is growing evidence that environmental 
and genomic stressors may induce significant 
levels of developmental instability (DI) 
(Palmer and Strobeck 1986; Palmer 1996; 
M0ller and Swaddle 1997; Pertoldi et al. 
2006a; Kurbalija et al. 2010), suggesting that 
DI could serve as an early signal in 
monitoring the effects of stress in populations 
(Leary and Allendorf 1989). 
Developmental stability is the ability of an 
organism to buffer environmental and genetic 
disturbances which affect the developmental 
capacity of a particular phenotype (Zakharov 
1981; Palmer 1996). The failure to correct for 
random accidents during development may be 
manifested as fluctuating asymmetry (FA), 
which is defined by small deviations from 
perfect bilateral symmetry of morphological 
structures (Leary and Allendorf 1989; Palmer 
and Strobeck 1986). 

The genetic basis of developmental stability 
has been much debated over the last five 
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decades. The increase or decrease of DI as a 
consequence of the genomic stress has been 
explained by two major hypotheses. The first 
argues that the level of stability is a reflection 
of the underlying genomic heterozygosity 
(Lerner 1954), while the second argues that 
stability reflects the general level of genomic 
coadaptation (Dobzhansky 1950). 

The coadaptive gene complexes are 
established over the evolutionary history of 
the genome via natural selection and are 
defined as a specific balance between loci 
within the genome (Dobzhansky 1950). Such 
coordination within the genome protects 
individuals from developmental accidents, 
which can be caused by both environmental 
and genetic factors (Parsons 1990). A 
breakdown of coadaptation can be manifested 
in an individual as a decreased ability to 
develop an optimal phenotype due to 
increased DI (Leary et al. 1985; Markow 
1995). 

Heterozygosity theory predicts that the level 
of heterozygosity is correlated with DI 
(Lerner 1954; Livshits and Kobyliansky 
1985), such that heterozygous individuals are 
better adapted to genomic and environmental 
perturbations than homozygous individuals 
due to higher developmental stability. The 
theory also states that single loci coding for 
enzymes related to metabolic efficiency 
influence developmental stability of different 
morphometric characters (Palmer 1996; 
M0ller and Swadle 1997), predicting that 
levels of heterozygosity at loci coding for 
functional proteins will be inversely correlated 
with level of DI (Lerner 1954; Livshits and 
Kobyliansky 1985). However, an important 
assumption is that heterozygosity of various 
genetic markers accurately reflects the 
heterozygosity of the entire genome, or at 
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least the heterozygosity of the loci that 
contribute to the formation of the 
morphological phenotype (Mitton and Grant 
1984; Chakraborty 1987). 

Whether heterozygosity and/or genomic 
coadaptation influence parameters of 
developmental stability is still unclear 
(Pertoldi et al. 2006b). Available data indicate 
a tendency of FA to increase with inbreeding 
(Palmer and Strobeck 1986; Waldmann 1999, 
2001; Babbitt 2006; Fessehaye et al. 2007), 
hybridization between species (Hochwender 
and Fritz 1999; Stamenkovic-Radak et al. 
2009), or between populations (Waldmann 
1999; Kurbalija et al. 2010), although several 
studies report exceptions to this pattern 
(Clarke et al. 1992; Sheridan and 
Pomiankowski 1997; Trotta et al. 2005). A 
consensus does appear to exist in evidence 
that any stress would increase phenotypic 
variance (V p ) of most quantitative traits 
(Bijlsma and Loeschcke 1997; Hoffmann and 
Parsons 1997; Hoffmann and Hercus 2000; 
Pertoldi et al. 2001a,b). 

It is well known among Drosophila spp. that 
V p of individuals collected in the field is 
greater than those reared in a laboratory 
setting (Coyne and Beecham 1987; Imasheva 
et al. 1994; David et al. 1997; Gibert et al. 
1998; Gibert et al. 2004). The higher 
phenotypic variability in flies from the wild is 
generally attributed to higher environmental 
variance, resulting in lower heritability 
(Coyne and Beecham 1987; Gibert et al. 1998; 
Hoffmann and Hercus 2000). It has been 
suggested that several mechanisms can alter 
the components of V p (Waddington 1942), 
and recent studies have aimed to identify their 
respective roles within an evolutionary 
context (e.g. Zakharov 1992; Debat and David 
2001; Meiklejohn and Hartl 2002). 
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Sternopleural bristle number is a widely 
studied characteristic in quantitative genetics 
(Mackay 2004) that displays a rapid response 
to artificial selection in Drosophila spp. 
(Ramuson 1955; Barnes and Kearsey 1970; 
Schnee and Thompson 1984). The genetic 
variance of SB obtained from IF analysis in 
Drosophila melanogaster was higher at 
extreme temperatures (Imasheva et al. 1998; 
Petavy et al. 2006). Some studies show a 
positive correlation between FA of 
sternopleural bristle number and temperature 
variations (Parsons 1961; Imasheva et al. 
1997; Bubliy et al. 2000; Petavy et al. 2006) 
and nutritional variations as stress conditions 
(Imasheva et al. 1999). The results of these 
studies support the rather popular point of 
view that environmental stress is associated 
with higher FA levels (M0ller and Swaddle 
1997). On the contrary, Woods et al. (1999) 
reported that stress did not induce an increase 
in the V p or FA of sternopleural bristle 
number in D. melanogaster. Despite 
heterogeneity among results in the association 
between different kinds of environmental 
stresses and FA and V p as estimators of DI, 
evidence is scarce regarding the effects of 
interspecies or interpopulation hybridization 
as a genomic stressor and its correlation with 
FA of this meristic trait in Drosophila species 
(Gilligan et al. 2000; Andersen et al. 2002; 
Vishalakshi and Singh 2008). 

Inversion polymorphism of different 
Drosophila species was used as a model 
system for studying processes involved in 
adaptation and genetic diversity (Balanya et 
al. 2004; Hoffmann et al. 2004; Stamenkovic- 
Radak et al. 2008). As crossing-over is 
suppressed within the inversion loops of 
heterokaryotypes, all genes within the 
inverted segments segregate as a linked group 
representing one physical and functional unit, 
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called the supergene, the different 
arrangements of which can be referred to as 
allelic complexes (Krimbas 1993; Rasic et al. 
2008; Santos 2009). D. subobscura is a 
Palearctic species that displays rich inversion 
polymorphism in all five acrocentric 
chromosomes of the set (Krimbas and Loukas 
1980; Krimbas 1992, 1993). Assuming a 
relatively long-time of selection on the linked 
genes within inverted regions, D. subobscura 
represents a suitable species for testing the 
heteozygosity vs. coadaptation hypotheses. 

This paper focuses on the coadaptive aspect of 
inversion polymorphism in D. subobscura 
populations from three ecologically and 
topologically distinct habitats, which possess 
a certain degree of genetic differentiation due 
to their different evolutionary histories. The 
aim of the study is to investigate the level of 
DI estimated by FA and V p in sternopleural 
bristle (SB) number in intra- and 
interpopulation hybrids through 3 generations 
(P, Fl and F2) in isofemale lines (IF) of D. 
subobsura. The study intends to contribute to 
understanding the effects of coadaptations and 
heterozygosity on the level of developmental 
homeostasis in intra- and interpopulation 
hybrids. 

Materials and Methods 

Population samples 

For the present study, D. subobscura 
individuals were sampled in Serbia using 
fermented fruit traps. The flies were collected 
from three localities (beech-B, oak-0 and 
Botanical Garden-BG). The first locality, 
Beech wood (B) (Abieto-fagetum), is situated 
between N 43°33' 28.43" and E 20°45' 10.96" 
(Mountain Goc, Central Serbia). The second 
is Oak wood (O) (Fraxineto-quercetum) 
situated between N 43°32' 57.38" and E 
20°40' 2.32" (Mountain Goc, Central Serbia). 



Kurbalija Novicic et al. 

These two woods have distinctive 
microclimates; Beech wood has the highest 
altitude (875 m) and the highest humidity with 
great vegetation coverage, while Oak wood 
(787 m) has more sparsely distributed trees 
and is slightly warmer. The third locality is 
the Botanical Garden at the University of 
Belgrade (BG) (Arboretum- Corilus colurna, 
Celtis australis), situated between N 44° 49' 
and E 20 ° 28' at 87 m above sea level, 
representing a more urbanized environment 
with a unique microclimate and high 
anthropogenic influence. 

The females collected at these three localities 
were individually used to obtain isofemale 
lines (IF) reared on the common cornmeal- 
sugar-yeast-agar medium for Drosophila. The 
progeny of these IF lines were used as the 
parental generation in the experiment in order 
to preserve the high amount of genetic 
variability from natural population. All 
cultures were maintained and all experiments 
performed under constant laboratory 
conditions at 19°C, approximately 60% 
relative humidity, light of 300 lux and 12:12 
L:D cycles. 

Population genetic structure and data 
analysis 

Analysis of inversion polymorphism was 
carried out for the wild captured D. 
subobscura males, which were individually 
crossed with virgin females from Kiisnacht 
laboratory stock, homozygous for Standard 
gene arrangement at all five large 
chromosomes. Salivary glands from third- 
instar larvae were squashed and chromosomes 
stained with aceto-orcein solution. Eight 
larvae were analyzed from the progeny of 
each of the crosses performed. The 
chromosome map from Kunze-Muhl and 
Muller (1958) was used for the cytological 
analysis of gene arrangements. The 
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designation of gene arrangements followed 
that of Kunze-Muhl and Sperlich (1955). The 
analysis included in total 56 males (112 
autosomes, 56 sex chromosomes) from Oak 
population, 44 males from Beech population 
(88 autosomes, 44 sex chromosomes), and 52 
males from Botanical Garden (104 autosomes, 
52 sex chromosomes). 

Z-statistics (Zar 1999) were used to assess the 
differences between frequencies of gene 
arrangements in the pairs of analyzed 
populations. The G-test (Sokal and Rohlf 
1980) was used to determine population 
subdivision by determining the homogeneity 
of gene arrangement frequencies between 
pairs of populations from different forest 
communities on all chromosomes and 
autosomes. 

Experimental design 

The experiment used 63 IF lines from Oak, 38 
from Beech, and 64 from the Botanical 
Garden; the progeny of these IF lines were 
used as the parental (P) generation in the 
experiment. Virgin males and females were 
separated within each IF line upon emerging 
and intra- and interpopulation crosses were 
made four days after eclosion. 

Intrapopulation (B x B, O x O, BG x BG) and 
interpopulation crosses (B x O, BG x O, BG x 
B) were made among IF lines of the three D. 
subobscura populations. Both direct and 
reciprocal crosses were made in order to take 
into account any maternal effect (Table 1). 
The progeny (six males, six females) from 
each cross were transferred to fresh vials to 
obtain Fl and F2 generations. Individuals in 
P, Fl, and F2 generations from intra- and 
interpopulation crosses (B x B, O x O, BG x 
BG, B x O, BG x O, BG x B) were frozen at 
-20°C and used for further sternopleural 
bristles counting procedure. 
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Count of sternopleural bristle number 

The left and right side of the sternopleural 
region from each fly was observed under a 
100x binocular microscope and the SB 
number was scored and counted twice, first on 
the right side (R), then on the left side (L). 
This meristic characteristics and others may 
be measured without error thus avoiding the 
need for replicate measurements (Palmer 
1994). According to most other similar 
investigations, all short, medium-sized, and 
long bristles were counted. 

Fluctuating asymmetry data analyses 

Before interpreting FA population estimates, 
several statistical assumptions were made. In 
the folowing statistical analysis, males and 
females from each population, each 
generation, and direct and reciprocal crosses 
were analyzed separately. Because of the large 
number of tests conducted, the sequential 
Bonferroni test (Rice 1989) was used. 
Fluctuating asymmetry is characterized by 
normal distribution of the right-side minus 
left- side differences with a mean of zero and a 
normal distribution (Palmer and Strobeck 
1986). No deviations from normal distribution 
was found using Shapiro-Wilk (W) and Chi- 
squared (% ) tests.. One-sample J-tests were 
performed to test for a departure from the 
mean of (R - L) from the expected mean of 
zero; no directional asymmetry was detected 
in all tested samples. Tests for correlation 
between SB number and FA were done 
because significant correlation may affect 
results and interpretation (Palmer 1994; 
Palmer and Strobeck 2003). Linear 
dependence of FA on the mean SB number 
was tested by linear regression; FA acted as 
the explanatory variable and SB number as the 
dependent variable. In 43.05% cases we found 
significant positive corelation between FA and 
SB number (Table 2). 
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Because higher variance was typically 
observed in cases where more bristles were 
found, a log transformation was used on 
variance values and further tests were done 
with transformed values to account for a 
scaling effect of the mean. The FA1 index 
(Palmer 1994) was measured as IR - LI 
between sides in all samples of 
intrapopulation and interpopulation hybrids, 
both direct and reciprocal crosses, separately 
for males and females through P, Fl, and F2 
generations. The FA1 index is the one of the 
most frequently used indices for describing 
the level of FA in a sample. It also acts as an 
unbiased estimator of the sample standard 
deviation, and is recommended for testing FA 
differences between three or more samples 
(Palmer and Strobeck 1992). However, 
because of the significant presence of positive 
correlation between FA and SB number, we 
also used FA2 index (Palmer 1994). The FA2 
index (mean (IR - LI/((R+L)/2))) is commonly 
used only where clear evidence exists of a size 
dependence of IR - LI among individuals 
within a sample. 

The multivariate ANCOVA was conducted to 
determine effects of intra- and interpopulation 
hybridization, generation, sex as an 
independent variable, and direct and 
reciprocal crosses as covariate variables while 
controlling for FA1, FA2, and mean 
sternopleural bristle number, the latter acting 
as a dependent variable. The analysis was 
performed using general multivariate linear 
model in SPSS software. 

Student J-tests and F-tests were used to test 
differences in mean and variances of the 
sternopleural bristle number and FA between 
sexes, populations, generations, and type of 
cross. All tests treated sex, population, 
generation, and type of cross as separate 
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variables. The conservative F-test and t-test 
were used to reduce the possibility of a Type 1 
error. Furthermore, we have performed 
corrections for multiple comparisons and 
therefore our results are both conservative and 
robust, and not affected by eventual 
interactions between the different factors. All 
statistical analyses were performed using 
PAST software (Hammer et al. 2001). 

Results 

Population genetic structure 

Frequencies of gene arrangements and 
parameters of inversion polymorphism 
(Degree of Heterozygosity, Inversion Density, 
Index of Free Recombination) on five 
acrocentric chromosomes of D. subobscura 
population from three different natural 
habitats are shown in Table 3. Inversion 
polymorphism analysis of Beech population 
(B) detected 15 stuctural chromosomal types 
and 13 inversions; in the Oak (O) and the 
Botanical garden populations (BG), 16 
stuctural chromosomal types and 15 
inversions were detected. 

Results of the G-test gave the level of 
interpopulation differences in gene 
arrangement distribution for each individual 
chromosome. The results showed significant 
differences in gene arrangement distribution 
between B and O populations for four 
chromosomes: A (G = 10.91, p < 0.01), J (G = 
14.49, p < 0.01), U (G = 11.10, p< 0.01), and 
E (G = 9.52, p < 0.05), and for all 
chromosomes in total (G = 49.24, p < 0.01) 
(Table 4). 

Interpopulation differences in gene 
arrangement distributions of each 
chromosome were found between B and BG 
populations: U (G = 9.70, p < 0.01), E (G = 
12.84, p < 0.01), O (G = 14.44, p < 0.01), and 
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for all chromosomes in total (G = 39.25, p < 
0.01) (Table 4). 

Interpopulation differences in gene 
arrangement distributions on each 
chromosome were found between O and BG 
populations: J (G = 17.02, p < 0.01), U (G = 
30.44, p < 0.01), E (G = 31.23, p < 0.01), O 
(G = 24.53, p < 0.01), and for all 
chromosomes in total (G = 116.95, p < 0.01) 
(Table 4). 

The results of the Z-test showed the level of 
interpopulation differences in individual gene 
arrangement frequencies and revealed 
significant differences in some gene 
arrangement frequencies between each pair of 
tested populations (B and O, B and BG, O and 
BG) (Table 5). 

The results showed significant differences in 
some gene arrangement frequencies for all 
chromosomes of the set between some 
populations (B and O, BG and O). There were 
significant differences in gene arrangement 
frequencies between B and BG populations, 
for U, E and O chromosomes (Table 5). Gene 
arrangement Ei +2+ 9+i2 was detected only in O 
population, as well as gene arrangement 06 
which was detected only in BG population 
with a frequency of less than 3%. 

Multivariate testing of fluctuating 
asymmetry and mean of SB number 

The results of the multivariate test using a 
general linear model are presented in Table 6. 
Significant differences were detected between 
populations; intrapopulation and 

interpopulation hybrids (F = 29.524, p < 
0.01), generations (F = 113.114, p < 0.01), 
and type of cross (F = 2.995, p < 0.05). 
Significant interactions were also found 
between population x generation (F = 59.919, 
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p < 0.01) and generation x sex (F = 11.022, p 
<0.01). 

The results test of between-subject effect are 
shown in Table 7. The significant differences 
were found between populations (FA1, F = 
22.471, p < 0.01; FA2 ,F = 20.153, p < 0.01; 
mean SB number, F = 19.895, p < 0.01), 
generations (FA1, F = 49.513, p < 0.01; FA2, 
F = 85.854, p < 0.01; mean SB number, F = 
65.351, p < 0.01). Significant interaction was 
also found between population x generation 
(FA1, F = 15.901, p < 0.01; FA2, F = 13.197, 
p < 0.01; mean SB number, F = 44.140, p < 
0.01). 

Univariate testing of fluctuating asymmetry 
and mean of SB number 
Intrapopulation hybridization. 

Changes of mean and variance across 
generations in males. Analysis of the 
difference in mean and variance of the 
average SB number in males is given in Table 
8a. Generally, very low variability in 
sternopleural bristle number in D. subobscura 
populations was observed. A significant 
increase of the average SB number was 
detected through generations (P < Fl < F2) in 
OxO hybrids resulting from direct cross. In 
B x B hybrids, a significant increase in SB 
number was detected in the Fl generation 
when compared to the parental generation, 
followed by a decrease in SB number in the 
F2 generation in both direct and reciprocal 
crosses. 

In BG x BG hybrids, a trend of constant 
increase in SB number was detected across 
generations (P < Fl < F2) in both direct and 
reciprocal crosses. 

No significant change in variance across 
generations was detected except in the case of 
BG x BG hybrids for both direct (F P , F i = 1.58, 
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p < 0.05; Ffi, F2 = 1.45, p < 0.01) and 
reciprocal crosses (F PjF i = 1.69, p < 0.05; 
F F i, F2 = 1.32, p< 0.05). 

Changes of mean and variance across 
generations in females. Analysis of the 
difference in mean and variance of the 
average SB number in females is given in 
Table 8b. The statistically significant increase 
of SB number was successive through 
generations for all hybrid groups in both direct 
and reciprocal crosses. Regarding the 
variance, no statistically significant change 
was detected across generations in all tested 
groups of intrapopulation hybrids for females. 

Changes of FA across generations in males. 

Analysis of the FA1 index between 
generations for SB number in males shows no 
significant difference between generations in 
the direct cross. There was, however, 
significant decrease of FA value across 
generations in O * O (t P , F i = 3.48, p < 0.01; t 
P ,F2 = 4.00, p < 0.01) and B x B (t Fi>F2 = 2.41, 
p < 0.05) hybrids in reciprocal crosses (Table 
9a). 

Analysis of the FA2 index between 
generations for SB number in males shows no 
significant difference between generations in 
the direct cross, although there was significant 
decrease of FA value across generations in O 
x o (t P ,Fi = 2.84, p < 0.01; t P , F2 = 3.66, p < 
0.01) hybrids in reciprocal crosses (Table 8a). 

Changes of FA across generations in 
females. Analysis of FA differences between 
generations using FA1 index for SB number 
in females shows no significant difference 
between generations in both direct and 
reciprocal crosses except in the case of B x B 
hybrids (t P>F i = 2.03, p < 0.05) in direct cross 
and O x O hybrids (t FlfP2 = 3.38, p < 0.01) in 
reciprocal cross (Table 9b). 
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Analysis of FA differences between 
generations using FA2 index for SB number 
in females shows no significant difference 
between generations in both direct and 
reciprocal crosses, except in the case of O x O 
hybrids in reciprocal cross, where FA value 
significantly decreased across generations (t 
fi,f2 = 3.69, p < 0.01). In the case of B x B 
hybrids, a constant decrease of FA across 
generations was detected (t PjF 2 = 2.60, p < 
0.05; t fi,f2 = 2.41, p < 0.05) (Table 10b). 

Inter-population hybridization 

Changes of mean and variance across 
generations in males. Analysis of the 
difference in mean and variance of the 
average SB number in males is given in Table 
8a. The results showed significant differences 
in SB number and variance in most 
interpopulation hybrids groups in both direct 
and reciprocal crosses. In B x O hybrids an 
increase in SB number in Fl generation was 
detected when compared to the parental 
generation, followed by a significant decrease 
in the F2 generation (t fi,f2 = 2.35, p < 0.05) in 
direct crosses. A different trend was found for 
reciprocal crosses (P > Fl > F2). In the case of 
BG x O hybrids, the significant increase in 
average SB number across generations was 
detected (P < Fl < F2) in direct crosses. In BG 
x B hybrids, a significant decrease of SB 
number in the Fl generation was found when 
compared to the parental generation, followed 
by a significant increase in SB number in the 
F2 generation for both direct and reciprocal 
crosses. 

The variance in SB number in B x O hybrids 
showed a statistically significant increase in 
F2 compared to Fl generation in both direct 
(F fi,f2 = 1-83, p < 0.01) and reciprocal 
crosses (F fi,f2 = 1.16, p < 0.01). In the case of 
BG x O hybrids, the variance showed a 
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significant constant increase across 
generations (P < Fl < F2) for both direct and 
reciprocal crosses. In the BG x B hybrids a 
significant increase of variance was detected 
in F2 when compared to Fl generation in 
reciprocal crosses (F PjF2 = 1-82, p < 0.05; F 

F1,F2= 1.61, P < 0.01). 

Changes of mean and variance across 
generations in females. Analysis of the 
difference in mean and variance of the 
average SB number in females is given in 
Table 8b. In B x O hybrids an increase of 
average SB number was detected in the Fl 
generation when compared to the parental 
generation followed by a decrease in the F2 
generation in reciprocal crosses. In BG x O 
hybrids a statistically significant increase of 
SB number was detected across generations (P 

< Fl < F2) in both direct and reciprocal 
crosses. In BG x B hybrids a significant 
decrease in average SB number was detected 
in the Fl generation when compared to the 
parental generation, and an increase in the 
same value for the F2 generation (P > Fl < 
F2) in both direct and reciprocal crosses 

A decrease in variance was detected in the Fl 
generation followed by a significant increase 
in the F2 generation for B x O hybrids in both 
direct (F P , F2 = 1.71, p < 0.05; F Fi , F2 = 2.03, p 

< 0.01) and reciprocal crosses (F PjF 2 = 2.73, p 

< 0.01; F fi, F2 = 2.93, p < 0.01). In BG x O 
and BG x B hybrids a significant increase of 
variance was detected in the F2 generation 
when compared to the Fl generation in both 
direct and reciprocal crosses. 

Changes of FA across generations in males. 

Analyses of the FA1 index between 
generations for sternopleural bristles in males 
are given in Table 9a. The results showed a 
statistically significant decrease in FA values 
across generations in B x O hybrids in both 
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direct (t P , F2 = 4.50, p < 0.01; t F i >F2 = 6.03, p < 
0.01) and reciprocal crosses (t PF 2 = 4.80, p < 
0.01; t fi, F2 = 7.07, p < 0.01). In BG x O 
hybrids a significant decrease of FA was 
detected in the F2 generation when compared 
to the Fl generation in both direct (t fi,f2 = 
2.19, p < 0.05) and reciprocal crosses (t F i, F2 = 
2.86, p < 0.01). In BG x B hybrids a 
significant increase of FA was detected across 
generations in both direct (t PjF2 = 2.57, p < 
0.05) and reciprocal crosses (t PjF i = 2.44, p < 
0.05). 

Analyses of the FA2 index between 
generations for SB number in males are given 
in Table 10a. The results showed a 
statistically significant decrease of FA values 
across generations in B x O hybrids in both 
direct (t P , F2 = 5.41, p < 0.01; t F i, F2 = 6.44, p < 
0.01) and reciprocal crosses (t PjF 2= 4.67, p < 
0.01; t fi,f2 = 7.50, p < 0.01). In BG x O 
hybrids a decrease of FA values was detected 
in the F2 generation when compared to the Fl 
generation in reciprocal crosses (t F i jF 2 = 3.17, 
p < 0.05). In BG x B hybrids a decrease of FA 
values was detected across generations in 
reciprocal crosses (t P , F i = 2.22, p < 0.05; t PjF 2 
= 2.60, p< 0.01). 

Changes of FA across generations in 
females. The results of FA1 analyses between 
generations for SB number in females are 
given in Table 9b. The results showed a 
significant increase in FA for B x O hybrids in 
the F2 generation when compared to the Fl 
generation in both direct (t F i jF2 = 7.16, p < 
0.01) and reciprocal crosses (t F i iF2 = 6.59, p < 
0.01). In BG x O hybrids a significant 
decrease of FA was found in reciprocal 
crosses (t F i >F2 = 3.19, p < 0.01). In BG x B 
hybrids a significant decrease of FA was 
found across generations in direct crosses (t 
p, F2 = 2.57, p< 0.05). 
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The results of FA2 analyses between 
generations for SB number in females are 
given in Table 10b. The results showed an 
increase in FA value in the Fl generation 
when compared to the parental generation, 
along with a significant decrease in the F2 
generation for B x O hybrids (P < Fl > F2) in 
both direct and reciprocal crosses. In BG x O 
hybrids a significant decrease of FA was 
found in the F2 generation when compared to 
the Fl generation in direct crosses (t F i,F2 = 
3.08, p < 0.01), while a statistically significant 
decrease in FA across generations was found 
in reciprocal crosses (t fi,f2 = 2.63, p < 0.01 ;t 
fi,f2 = 3.68, p < 0.01). In BG x B hybrids a 
statistically significant decrease of FA was 
detected across generations in direct crosses (t 
P ,F2 = 3.10, p < 0.01; t fi,f2= 3.58, p < 0.01). 
The same trend of FA variability across 
generations was also found in reciprocal 
crosses (t F i,F2= 1.97, p < 0.001). 

Discussion 

The results showed that populations of D. 
subobscura from three ecologically and 
topologically distinct habitats possess a 
certain degree of genetic difference, likely due 
to their different evolutionary histories. The 
gene arrangements are carriers of various 
alleles that are differently favored in diverse 
environmental conditions and prove in most 
cases to be the major factor determining the 
gene arrangement frequencies in natural 
populations of D. subobscura (Andjelkovic et 
al. 2003). The results presented here further 
emphasize the relationship between gene 
arrangement frequencies and local adaptation 
to ecologically different microhabitats in D. 
subobscura. 

The inversion polymorphism analysis showed 
that analysed populations differ in the 
frequencies of some gene arrangements that 
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are in agreement with previous results 
(Andjelkovic et al. 2003; Stamenkovic-Radak 
et al. 2008; Jelic et al. 2009). In general, the 
pattern of inversion polymorphism of the 
three analyzed D. subobscura populations was 
consistent with the observed inversion 
polymorphism of D. subobscura in the 
southeast margin of the Central European 
region (Krimbas 1993). Therefore, the results 
suggest that the three populations in this study 
are suitable for testing the coadaptation vs. 
heterozygosity hypothesis. 

This paper focuses on the coadaptive aspect of 
genetic variability at the population level and 
its relationship to interpopulation 
hybridization as a genomic stressor. The trait 
used in the study was a typical meristic 
character widely considered as a threshold 
trait (Falconer and Mackay 1996; Mackay 
2001, 2004). Some data suggest the existence 
of particular genes that affect the genotype x 
environmental covariance of SB number 
(Kristensen et al. 2005; S0rensen et al. 2007). 
The underlying continuous variable is both 
genetic and environmental in origin and could 
be measured and studied equally as a metric 
character but with greater precision because of 
non-existence of measurement error. 

The results regarding the mean of SB number 
for intra- and interpopulation hybrids of D. 
subobscura cannot explicitly reveal the 
significance of either of the two hypotheses, 
since no general trend has emerged for either 
intra- and interpopulation hybrids. However, 
the observed variance of the SB number 
across generations gives a completely 
different insight. 

The variance of SB number after intra- and 
interpopulation hybridization showed that 
intrapopulation hybrids have no significant 
change of variance across generations in most 
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of the hybrid groups, compared to 
interpopulation hybrids. However, a 
significant change of variance across 
generations is observed in interpopulation 
hybrids. A reduction of V p in the Fl 
generation compared to the parental 
generation was obtained in most of the 
crosses, signifying that the hybrids had 
reached a higher level of homeostasis in the 
Fl generation, likely due to higher 
heterozygosity. Also, a significantly higher 
level of V p in the F2 generation was obtained 
when compared to the Fl generation, probably 
due to disruption of coadapted gene 
complexes, which can be attributed to 
crossing over in formation of most of the 
genotypes in the F2 generation. 

The analysis of the variability of FA of SB 
number using FA1 and FA2 indices (Palmer 
1994) after intra- and interpopulation 
hybridizations showed that intrapopulations 
hybrids are relatively stable. Analysis showed 
no significant differences in FA in most of the 
cases between generations when compared to 
the interpopulation hybrids across generations. 
This result was expected when the effects of 
intrapopulation crosses are taken into account, 
which simulate random mating in natural 
populations. These crosses represent control 
conditions for genomic stress. On the other 
hand a significant trend of change in FA was 
observed across generations in interpopulation 
hybrids, which was not observed in 
intrapopulation hybrids. 

Theoretically, a heterosis in the Fl generation 
and outbreeding depression caused by 
disruption of coadaptive gene complexes in 
the F2 generation for interpopulation hybrids 
was expected, but results showed a post- 
hybridization decrease in FA for 
interpopulation hybrids in both Fl and F2 
generations. 
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The results showed a significant decrease of 
FA through generations each with the lowest 
FA detected in the F2 generation in all 
interpopulation hybrids. Therefore, it is likely 
that the observed difference between intra- 
and interpopulation hybrids regarding level of 
FA variability is due to increased 
heterozygosity in the Fl generation in 
interpopulation hybrids, and thus the 
increased developmental homeostasis and 
reduced FA. An increase of heterozygosity 
also reduced the FA in the F2 generation, 
which could be a result of selection pressure 
acting on F2 genotypes. Selection could 
eliminate developmentally unstable and 
extreme phenotypes from the population at 
very early stage and overall hybrid fitness 
could increase (Polak et al. 2002). The 
increase of FA in the F2 generation compared 
to the Fl generation was not obtained, 
indicating the possibility that disruption of 
coadaptation occurred and the less stable 
genotypes produced by crossing over were 
selected against. Therefore, in the F2 
generation, we obtained only individuals with 
suitable gene combinations and higher 
developmental stability. This suggests that the 
most probable mechanism providing 
developmental homeostasis in a population is 
the heterozygote or hy bride superiority, also 
known as overdominance. 

Heterosis can be considered the opposite of 
inbreeding depression; both mechanisms share 
the same underlying causes (Crow 1993). The 
primary cause of heterosis is the sheltering of 
deleterious recessive alleles in hybrids. In 
addition, increased heterozygosity would 
increase the fitness of hybrid individuals at 
loci where heterozygotes have a selective 
advantage over homozygote types. A 
subdivision of natural populations can provide 
an appropriate condition for heterosis, as 
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different deleterious recessive alleles could 
drift to relatevely high frequencies in different 
subpopulations. Therefore interpopulation 
hybrids produced by mating between 
immigrant individuals are expected to have 
greater fitness than residents (Whitlock et al. 
2000; Morgan 2002). 

The evidence for heterosis without 
coadaptation in Drosophila ananassae was 
presented in Sing (1972, 1985), showing that 
interracial hybridization does not lead to 
breakdown of heterosis associated with 
cosmopolitan inversions. These results 
provide observations that may carry weight on 
testing heterozygosity vs. coadaptation , but 
the findings of heterosis in interpopulation 
hybrids does not exclude hybrid breakdown in 
subsequent generations (Lynch and Walsh 
1997). 

In response to genomic stress, no significant 
difference was found between the majority of 
direct and reciprocal crosses, suggesting that 
parental population origin had no significance 
and that both sexes equally contribute to the 
inheritance of SB number. This suggests that 
sternopleural bristles, as a bilaterally 
symmetrical morphological trait, might be a 
good and consistent characteristic for 
measuring the level of developmental stability 
under genomic stress without a confounding 
maternal effect. These results contrast with 
wing FA in a previous study (Kurbalija et al. 
2010). 

There may be several explanations for the 
observed difference in the response of 
different bilateral symmetrical traits to 
genomic stress conditions. The development 
of distinct traits is probably under control of 
different gene complexes (Parsons 1990). A 
stress factor may also be specific for particular 
metabolic pathways and may not affect the FA 
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of all traits in the same way (Palmer and 
Strobeck 1986). Furthermore, different traits 
are under different natural selection pressure 
that depends on functional importance of a 
given trait. Sternopleural bristles are not 
generally regarded as a fitness -related trait 
and we do not expect them to be a highly 
canalized trait. 

Despite the fact that different outcomes can be 
expected after interpopulation hybridization, 
depending on evolutionary histories of 
populations and trait specific response, our 
results confirmed that DI could be a useful 
indicator of genomic stress in populations. 
Further knowledge on the consequences of 
interpopulation hybridization between wild 
populations as well as the capacity of small 
populations to adapt to local environmental 
conditions is urgently needed. 

The integration of experimental, theoretical, 
and applied conservation genetics will 
contribute to understanding of the 
methodological and applied aspects of 
conservation genetics. It seems that 
Drosophila is a useful species to bridge the 
gap between theoretical computer simulation 
and studies on natural populations. 
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Table 2. Tests for linear dependence of FA on SB average number using 
linear regression, of intra- and inter-population crosses in P, Fl and F2 
generations in males and females (direct and reciprocal crosses). 







Males 


1 iiniiks 






direct 


reciprocal 


direct 


reciprocal 


Cms 


generation 


JV 


correlation 


r 


/V 


correlation 


r 


N 


rorrrlitlon 


r 


/V 


correlation 


P 


Inlra-popiilation 




























1' 


16 


•0.19 




16 


0.637 




16 


0.0 IX 




16 


0.124 




OxO 


H 


73 


11.367 




71 


-0.035 




62 


ii. r. 




67 


0.195 






F2 


173 


11.263 




132 


ii :r. 




199 


0234 




167 


0.117 






P 


1.1 


(1217 




13 


0533 




13 


0321 




13 


-0.2 
































BxB 


Fl 


61 


U .' 




58 


0.0X8 




62 


0.017 






0,045 






F2 


115 


11.044 




116 


0 112 




144 


0205 




I3X 


0,097 






1" 


49 


(1.073 




49 


0.169 




49 


0.175 




49 


0.126 




Ili.OK, 


Fl 


215 


0.147 




222 


0.167 




207 


0.309 




20X 


0.I8X 






F2 


593 


O.IK 




615 


0.099 




712 


0, 1 76 




6X6 


0,217 




Intrr-populatlon 




























r 


42 


-11 21 




42 


0.276 




42 


0.134 




42 


0.293 




Bx<> 


Fl 


195 


-0.005 




1X6 


0.051 




192 


•0.045 




1X0 


0,021 






i : 


433 


0.361 




410 


0.432 




5IX 


0.311 




493 


0.2X 






p 


46 


-0.053 




46 


0,106 




46 


0.299 




46 


ii |4< 
































B(ixO 


1 1 


215 


0.22 




203 


0.004 




198 


0.107 




18- 


0.2 






n 


444 


0.167 




372 


0.25X 




528 


0,289 




4X4 


0,233 






p 


44 


0.132 




44 


0.23 




44 


0.421 




44 


0.344 




IU.U1 


Fl 




(1.092 




215 


0.23 




202 


0.096 




190 


0,204 






F2 


5x2 


II DM 




5X2 


0 1X9 




673 


0.24X 




MX 


0 204 





O, Oak population; B, Beech population; BG, Botanical garden population. 
p<0.05 *, p<0.0l**, p<0.00l*** 



Table I . Number and type of intra- and inter-population crosses 
for direct and reciprocal crosses. 





Direct cross | Reciprocal cross 




Type of cross 


Inter- isofemale line crosses 


Number of crosses 


Infra-population 








OxO 


O; x O ' 


O ' X O : 


63 


BxB 


Bi x B 


B ' X B; 


71 


BGxBG 


BG . x BG ' 


BG ' x BG , 


8« 


Inter-population 








BxO 


B . x O ' 


B ill; 


82 


BGxO 


BG ; x O ■' 


BG ' x O ; 


81 










BCixB 


BG xB 


B(i ' X B ; 


73 



O, Oak population; B, Beech population; BG, Botanical Garden 
population 



Table 3. Gene arrangement frequency (%) and inversion 
polymorphism parameters of Drosophila subobscura populations 
from three different habitats. 



Gene 
arrangement 


O (N=56) 




B(N=44) 


BG (N=57.) 


n 


% 


n 


% 


n 


% 


\-i 


s -i 


39.29° o 


29 


65.91",, 


31 


59.62",, 


\ 


32 


57.14% 


II 


25.00% 


1" 


36.54% 


A. 


2 


3.57% 


4 


9.09' ,. 


2 


3.X5",, 


.1 


3 


2.68",, 


In 


IS. IS",, 


20 


19 25",, 


J, 


109 


97.32° u 


- : 


M s.2",. 


84 


SO. 77",, 
















U„ 


8 


7.14% 


9 


10.23% 


24 


23.08% 


u,„ 


25 


22.32° o 


37 


42.05°i 


50 


48.08% 


U : 


79 


70.54% 


42 


4^ 73% 


30 


28.85% 


E„ 


18 


16.07% 


15 


17.05% 


40 


39.42",, 


E» 


17 


15.18% 


26 


29.55% 


30 


28.85% 


E 


74 


66.07" o 


4" 


53.41% 


34 


31.73% 


E|. 


3 


:.<>v ,. 


n 


0.00",, 


11 


0.00' 


o 


12 


10.71°,, 


•l 


10.23% 


24 


24 04",. 


O, 


0 


11.00",, 


II 


o.oo",, 




0.9f,"„ 


O, , 


44 


39.29° o 


45 


51.14% 


56 


52.88% 


















47 


41.96% 


30 


34.09° o 


16 


15.38% 


0,. 4 . : 


9 


N.04".. 


4 


4.55% 


6 


6.73% 


HZ 


1.64 




2.48 




2.17 




1FR 


86.08 




80.29 




80.96 




ID 


7.27 




6.89 




5.88 





O = Oak population, B = Beech population, BG =Botanical garden 
population 

HZ, Degree of Heterozigosity; IFR, Index of Free Recombination; 
ID, Inversion Density 

N = number of gene arrangement; N = number of analyzed males 



ournal of Insect Science | www.insectscience.org 



18 



Journal of Insect Science.Vol. I I | Article I 13 



Kurbalija Novicic et al. 



Table 4. Chromosomal differences in arrangement frequencies 
between three analyzed habitats of Drosophila subobscura 
populations. 



Chromosomes 


B/O 


B/BG 


O/BG 


A 


10.91 ** 






J 


14,49 *** 


0.03 


17.02 *** 


U 


1 1.01 ** 


9.70 ** 


30 4 ] *** 


1. 


9.52 * 


12.84 ** 


31.23 *** 


o 




1 4 44 * * 


24.53 *** 


all 


49 24 *** 


39 ^5 *** 


116.95 *** 



O = Oak, B = Beech ,BG = Botanical garden. p<0.05 * p<0.0l **, 
p<0.00l *** 

G-test values are given only for significant comparisons. 



Table 5. Differences in gene arrangement frequencies between 
three analyzed habitats of Drosophila subobscura populations. 



Gene 
arrangement 


B/O 


B/BG 


O/BG 


A„ 


3.76 *** 




(-)2.99 ** 


A, 


(_)4 S9 *** 


3.03 ** 


A; 








J-l 


5.29 *** 




(-)5.58 *** 


J, 


(-)5.29 *** 


5.58 *** 


u„ 




(-)3.33 *** 


(-)4.66 *** 


u,. 2 


4.26 *** 




(-)5.62 *** 


u„ 2 . 6 


(-)4.66 *** 


3. 82 *** 


8.66 *** 


E« 




(-1-4.64 *** 


(-15.25 *** 


1-. 


44 *** 




(.)3.44 *** 


E| 1 ; 


(-)2.59 ** 


4 I 1 *** 


6.94 *** 


E, 


(-)2.20* 




2.38 * 


o„ 




(.)3.34 *** 


(-)3.45 


o„ 






(-)2.09 * 




2.3S * 




(-)3.03 ** 


(W, 




4 29 *** 


6.08 *** 


0, . 









O = Oak, B = Beech ,BG = Botanical garden. p<0.05 *, p<0.0l **, 
p<0.00l *** 

Z-test values are given only for significant comparisons, 
v / 



Table 6. Results of multivariate test using General linear model. 



Effect 


\ aluc 


F 


Hypothesis 
df 


F.rrurdf 


SlCTifkance 


Intercept 


Pill.n's Trace 


0 ".I 


1 i:s.i 


3 


14765 


0 


tt K ! amM.1 


0.269 


1 ',',59 9-0a 


3 


re 


0 


Hotclling's Traee 


2.715 


1 ns9.9-9j 


3 


14765 


0 


Roy's Largest Rool 


2.715 


13359.979a 


3 


14765 


0 


type ol LTOSS 


I'illai's Trace 


(MM 


2.995a 


3 


14765 


0.0.3 


Wilks' 1 ambda 


0.999 


2.995a 


3 


14765 


0.03 


Holelling's Trace 


(IIHII 


2.995a 


•■ 


14765 


0113 


Roy's Largest R,hii 


0.001 


: 90s., 


3 


14765 


III! 


population 


I'illai's Trace 


0.016 


16.268 


15 


44-.ll 


0 


Wilks' Lambda 


0.984 


16.295 


15 


40760,07 


0 


Holelling's Trace 


0.017 


16.314 


15 


44291 


0 


Roy 's ! argesl Rool 


II III 


29 524b 


5 


14767 


(1 


generation 


I'illai's Irate 


0.02.1 


57 


6 


29532 


0 


Wilks' Lambda 


0.977 


57 764a 


6 


29530 


II 


Holelling's Trace 


0.024 


5)1.064 


6 


2952(1 


0 


Roy 's I argest Riml 


0.023 


113.114b 


3 


14766 


0 


sex 


Pillai's Trace 


n 


1.636a 


'. 


14765 


11 179 


Wilks' Lambda 


1 


1.636a 


3 


14765 


0.179 


Holelling's Trace 


0 


1 636a 


l 


1 |-I,s 


0.179 


Roy's Largest Root 


0 


1.636a 


3 


14765 


0.179 


population X 
generation 


Pillai's Trace 


0.038 


19.134 


ill 


44KII 


0 


Wilks' Lambda 


0.962 


19.33 


-,u 


4»»S s~ 


0 


Holelling's Trace 


Mill 


l'i 522 


30 


44291 


0 


Roy 's Largest Root 


0.0.57 


53.919b 


Hi 


14767 


0 


population x 

sex 


Pillai's Trace 


ami 


1.03 


15 


44301 


0.42 


Wilks' Lambda 


0.999 


1.03 


15 


40760,07 


0.42 


Hotellme's Trace 


0.001 


1.03 


15 


44291 


0.42 


Roy 's Largest Root 


o.ool 


2 IS3b 


5 


i 1767 


11.113 


generation X 

sex 


generation x sex 


0.003 


6.904 


6 


29532 


0 


Wilks' Lambda 


0.997 


6.905a 


6 


29530 


0 


Holelling's Trace 


0.003 


6.906 


6 


29528 


0 


Roy 's Largest Root 


0.002 


11.022b 


3 


14766 


0 


population \ 
generation x 
sex 


Pillai's Trace 


0.002 


11 93o 


3(1 


44:u 


11 566 


Wilks' Lambda 




11.930 


30 


4333s 57 


11 s,„. 


llolcllmi-.'s Tract 


(HHI2 


11 910 


30 


44291 


11 566 


Roy's Largest Root 


0.001 


. 434b 


HI 


476- 


11 15* 



a Exact statistic; b The statistic is an upper bound on F that yields a 
lower bound on the significance level; c Design 

v J 



Table 7. Results of tests of between-subjects effects using 
genelar linear model (multivariate). 



Source 


Dependent 
Variable 


Type III 
Sum of 
Squares 


df 


Mean Square 


F 


Significance 


Corrected Model 


1 A 1 


245 SS3' 


v. 


6. S3 


13.941 


0 


1 \: 


(.61,1 




1, 1 5. 


15.323 


11 


1 k 1 1 : 


I'i'.o ;:i 


;., 


53.787 


-In 1159 


(i 


Intercept 


i \i 


5 1,, s.5.1 


1 


546.859 


1 1 1 (. ! 92 


0 


FA2 


14.441 


1 


14.441 


1 195.89 


0 


(R^L|2 


45229.161 


1 


45229.161 


38756.33 


0 


ty pe of cross 


FA 1 


0.049 


1 


0.049 


0.101 


0.751 


LA 2 


0.007 


1 


0.1X17 


0.574 


0.449 


(R+LV2 


9 373 


1 


9.373 


8.032 


0.005 


population 


LAI 


55.047 


5 


11.009 


22.471 


0 


FA2 


1.217 


5 


0.243 


20.153 


0 


IR 1.12 


116.086 


5 


23.217 


19.895 


0 


generation 


1 \ : 


45 5 1,, 


2 


24 2 5s 


•l'i 5, | 


(l 


LA 2 


2.073 




1.037 


85,854 


0 


(RI.I2 


152.532 




76.266 


65.351 


0 


sex 


1 A 1 


0.002 


1 


11.002 


(1003 


0.95s 


IA2 


0.003 


1 


11003 


0.219 


1164 


(HI 1 2 


5 ti)3 


1 


5 so'. 


1 62 


0032 


population x 
generation 


1 A 1 


" 905 


III 


7.79 


15.901 


0 


LA 2 


1 5.14 


HI 


11 is., 


1 |97 


11 


(R'L|2 


515.117 


HI 


51.512 


44.14 


0 


population x sex 


LAI 


4.91)2 


5 


0.998 


2.038 


0.07 


LA 2 


0.11 


5 


0.022 


1.828 


0.104 


(RI.I2 


2 IIS 


5 


11424 


11 363 


0,874 


generation x sex 


FA 1 


4.072 




: 036 


4.155 


(Hill. 


1 \: 


1111" 




0.035 


2. 882 


0 056 


(R'L|2 


35.P9 


5 


17.589 


15.072 


0 


population x 
generation x sex 


FA I 


4,64 


10 


0.464 


0.947 


0.488 


I \? 


11 1 17 


Hi 


11(114 


1 1 16 


II 11 


(R 1 1 2 


i ; sin 


in 


i (.5 


1 151 


(1 297 


Lrror 


FA 1 


7234.S45 


14767 


0.49 






FA 2 


178.322 


14767 


0.012 






(R+D2 


17233.29 


14767 


1.167 






Total 


FA1 


1.3775 


14804 








FA 2 


,. ■; 


14804 








(R*L)2 


623851.75 


14804 








Corrected T otal 


FAI 


74S0.72S 


14803 








FA 2 


is i .is ; 


14803 








l K ■ 1 


[.in." 1.1 1 


1451,: 









a R Squared = .033 (Adjusted R Squared = .03 I); b R Squared = 
.036 (Adjusted R Squared = .034); = R Squared = .101 (Adjusted R 
Squared = .099) 
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Table 8a. The mean and variance (log transformed) for sternopleural bristle number in males of intra- and interpopulation crosses across 
generations (direct and reciprocal crosses). 





\ | direct 




reciprocal 


cross | generation 


mean±SE 


Itesl 


P 


variance 


Ftest 


P 


N 


mean±SE 


ttest 


P 


variance 


Ftest 


P 


Intra-population 
































P 


16 


5.84±0.19 


P<Fl.t=-2.06 


• 


0.0039 


P>F1.F=1.25 




16 


6.50x0.27 


P>FI,t=1.29 




0.004 


P>F1,F=I86 




OxO 


Fl 


73 


(v 3D - ().()•■> 


P<F2.t=-2.17 


• 


0.0031 


P<F2.F=1.46 




71 


6.20±0.09 


P>F2.t=0.97 




ii.oo '-2 


P>F2,F=1.24 






F2 


173 


6.48 0,09 


FKF2.t=-1.19 




0.0047 


FKP2. F=1.83 




132 


6.29±0.07 


FKF2.t=-0.81 


0.0027 


F1>F2.F=1.20 






P 


1.3 


5.15x0.17 


P<Fl.t=-4.94 


*** 


0.003 1 


P<F1.F=1.67 




13 


5.65x0.24 


P<Fl,t=-3.8l 


«* * 


0.0048 


P<F1.F=1.07 




BxB 


Fl 


60 


6.72 = 0.14 


P<F2.t=-4.45 


*** 


0.0051 


P<F2,F=1.I8 




58 


6.89x0.14 


P<F2.t=-2.27 


* 


0.0051 


P-F2.F=1.09 






F2 


115 


6.19x0.07 


Fl>F2.t=.3.58 


«•* 


0.0036 


F1>F2.F=0.I2 




1 Id 


6.36x0.10 


Fl>F2.t=3.05 


** 


0 0052 


FKF2.F=I.0.3 






1' 


49 


5.76-0,1 5 


P<Fl.t=-3.05 


- r 


49 


P>F1.F=I.58 


* 


49 


5.86±0.I4 


P<Fl.t=-2.36 


* 


0.0058 


P>F1.F=1.69 


* 


BGxBG 


Fl 


221 


6.I5±0.05 


P<F2.t--5.70 


** * 


0.003 1 


P -F2.p-I.09 




226 


6. 1 8x0.05 


P-F2.t--5.01 


**• 


0.0034 


P>F2,F=1.28 






F2 


593 


1,1,4 (1(14 


1 1 i :.( -6.38 


««» 


45 


i i r:.i i 45 


** 


(•15 


0.66.0 1(4 


FKF2.t=-6.02 


• ** 


0.0045 


II 1 2.1 1 32 


* 


Intcr-population 
































I' 


42 


5.77±0.I3 


P<F1.I=-1.39 




0 004 3 


P>F1,F=1.44 




42 


6.33x0.24 


P>Fl.t=l.58 




0.00.37 


P<FI.F=I.I9 




B\( i 


Fl 


I'M 


5.95:0.05 


P -F2.I 0.11 




0.003 


P<F2.F=1.27 




190 


6.09*0.06 


P>F2.t=4.31 


*** 


(44 


P<F2,F=I39 






F2 


433 


5.76x0.05 


Fl-F2.t-2.35 


* 


0.0055 


FKF2.F-I.83 




409 


5.67x0.05 


Fl>F2.t=5.ll 


«** 


0.0051 


FKF2.F=I.I6 






P 


46 


5.97x0. 10 


P- Fl.t -3.05 


** 


0.003 1 


P-=F1.F=1.21 




46 


6.06x0.10 


P<Fl.t=-2.21 




0.0023 


P-=F1.F=1.94 




HG\<) 


1 1 


:i i 


6.40:0.11(1 


P- F2.t— 2.50 


* 


0.(103" 


P 1 2.1 2.1 1 


• 


203 


6.40x0.07 


P-F 1.1— 0.95 




0.0044 


P-F2.F-2.42 


**♦ 




F2 


44(. 


6.44x0.06 


FKF2.I=-0.38 




0.0066 


FKF2.F=I76 


*** 


373 


6.23x0.06 


F1>F2.1=1.72 




0.0055 


FKF2.F=I24 


• 




P 


44 


6.31 ±0.1 4 


P>F1. 1=3.41 




0.0037 


P<F1.F=1.17 




44 


6.39±0.I3 


P>Fl.t=0.99 




0.0034 


P<FI.F=1.I4 




BGxB 


Fl 


207 


5.82x0.06 


P<F2.t=-1.06 




0.0043 


P-F2.F-1.33 




215 


6.25±0.06 


P'F2.t=-l.l3 




0.0039 


P<F2.F-1.82 


* 




i : 


584 


6.49x0.05 


1 1 1 2.1 -7 "1 




0.0049 


FKF2.F=1.I4 




586 


6.61 0.05 


Fl- F2,t=-3.9I 


X** 


0.0063 


FKF2.F=I.61 





O = Oak, B = Beech ,BG = Botanical garden. p<0.05 * p<0.0l **, p<0.00l *** 

The F test was used for testing differences in variances and t test were used for testing differences in means of SB number. 



Table 8b. The mean and variance (log transformed) for sternopleural bristle number in females of intra- and interpopulation crosses across 
generations (direct and reciprocal crosses). 



cross | generation 






reciprocal 


Intra-population 


meaniSF; 


I lest 


P 


\ ariance 


F test 


P 


N 


meaniSK 


1 tesl 


P 


variance 


F lesl 


P 




P 


16 


6.19±0.16 


P>Fl.t=0.27 




0.002 


P<FI.F=I.62 




16 


6.03x0.22 


P<Fl.t--l.26 




0.0055 


P>FI.F=I72 




OxO 


Fl 


62 


6.13x0.10 


P 1 2.1 -1 10 




0.0032 


P 1 2.1 2.30 




67 


6.32 1) 1)9 


P 1-2.1 -2.47 


* 


0.0032 


PI2.I 157 






1 2 


200 


6.49±0.08 


II 1 2.1 -2.43 


f 


0.0046 


FKF2. 1 .42 




167 


6.6. 0.117 


Fl '-■-F2.t--2.4I 


* 


0.0035 


II I2.F 1.09 






P 


13 


5.6510.29 


P 1 l.i -2.40 


* 


0.0073 


P 1 I.I 1.09 




13 


5.46 0.17 


P l l.i -2.25 


* 


0.002 1 


P 11.1 1.87 




BxB 


1 1 


62 


6.5U0.15 


P<F2.t=-2.80 


*** 


0.0067 


P F2.F=I.42 




55 


6.58x0.13 


P = F2.t=-2.l6 


* 


0.0039 


P- F2.F-I.92 






F2 


145 


0.02:0. 10 


Fl I2.t -0 60 




0.(105 1 


Fl 1 2.1 1.30 




138 


6.55 0.08 


11 12.1 0 18 




0.0041 


II 12.1 1.0.3 






P 


49 


5.64:0.14 


P 1 1.1 -3.60 




0.0062 


P 1 I.I 1.61 




49 


5.05 11 1 3 


P Fl.t -1.24 




0.0051 


P 1 I.I 1,28 




BGxBG 


1 1 


208 


6. 1 610.06 


P<F2.tx-7.14 




0.0038 


P>F2.F=I.I5 




21.3 


6.12x0.06 


P<F2,t=-5.l8 




0.004 


P>F2.F=I.05 






1 2 


710 


6.90 ■ 0.04 


FKF2.t--8.21 




0.0054 


FKF2.F=I40 




685 


6.82x0.04 


FKF2,r=-8.12 




0.0019 


FKF2,= I.22 




Inter-population 
































P 


42 


6.0610.14 


P '1 1.1 o .34 




0.0043 


P I- I.I- 1.18 




42 


6.08 9! , 3 


P Fl.t -0.68 




0.0032 


P 1 1.1 1,07 




B\( I 


Fl 


192 


6.0110.60 


P -1-2.1-0.98 




0.0036 


P 1 2.1- 1.71 


* 


182 


6. 1 7x0.06 


P>F2.rl.47 




0.003 


P<F2.F'-2.73 






1 2 


5 Id 


5.8710.05 


Fl I2i 1.46 




0.0073 


FKF2.p-2.03 




49.3 


5.81x0.05 


Fl>F2.t-3.88 


*** 


0.0088 


111 2,1 2 9-; 


*** 




P 


46 


6.1910.10 


P<Fl.t=-0.59 




0.0025 


P 1 LI 1.51 




46 


6.02x0.10 


P<Fl.t=-2.27 


* 


0.0026 


P I- I.I 1.3.3 




BCixO 


Fl 


196 


6.28:0.06 


P<F2,t=-2.34 


• 


0.00.38 


P 1 2.1 2.28 


* 


185 


6.45x0.09 


P<F2.t=-2.94 


** 


0.00.3 5 


P 12.1 2.19 


** 




F2 


527 


6.61 0.05 


FKF2,t=-3.57 




0.0058 


FKF2,F=I.51 




482 


6.56 0.05 


FKF2.t=-1.00 


0.005.8 


11 12,1 1.64 


*** 




P 


44 


6 35 : 0.10 


P 1 1.1 2.30 


* 


0 6062 


P -F1.F=I .6I 


* 


4-1 


6.35x0.11 


P>Fl,t=1.49 




0.002" 


P<F1.F=I70 


* 


BGxB 


1 1 


207 


5.9810.07 


P-;F2.t-2.29 


* 


0.0038 


P>F2.F-I.I5 




191 


6.10x0.07 


P-F2.t--2.95 


** 


0.00-17 


P<F2.F-2.22 






F2 


669 


6.7810.05 


FKF2,t— 8.50 




0.0054 


FKF2.F-1.40 


*** 


648 


6.92x0.05 


FKF2.t--8.17 




0.0061 


FKF2.F-I.31 


■ 



O = Oak, B = Beech ,BG = Botanical garden. p<0.05 * p<0.0l **, p<0.00l *** 

The F test was used for testing differences in variances and t test were used for testing differences in means of SB number. 



Table 9a. The FA I index differences between generations and type of crosses for sternopleural bristles in males (direct and reciprocal 
crosses); t test were used for testing differences in mean of SB. 



FAI-IR-U 






Direct 


Reciprocal 


miss 


Beneratton 


N 


mianxSl 


1 lesl 


P 


N 


meanxSF 


1 lest 


p 


Intra-population 




















P 


16 


11.62 :0.22 


P-FI.1=0.2S 




16 


1 .12=0.22 


P 1 1,1 ' 4s 




OxO 


II 




1) 57 91 (17 


P 1 2.1 II 72 




71 


0.52x0.07 


P [ 2.1 4.00 






F2 


17! 


0.49x0.05 


11 -1-2.1 0.91 




132 


0.45x0.05 


II 1 2,1 -0.77 






P 


1.3 


0,46:0 IS 


I" 1 1.1 II 3" 




13 


0.69x0,13 


P III 'I 36 




IKH 


1 1 


60 


0 510(18 


p I :.i mis 




58 


0.62:0.09 


P 1 2.1 1 .94 






i : 


115 


II 45 -0 115 


1 1 1 2.1 0 85 




1 16 


040:005 


Fl -K2.I 2.41 






p 


49 


0.5 1 :0 II) 


P 1 1.1 -0 3| 




49 


0,45:0 12 


P 1 1,1 -1 72 




BOxBQ 


1 1 


221 


0.5410.03 


P K2.I -0.77 




226 


0.62x0.04 


P- 1 2.1 -1.49 






i : 


593 


0.5810.03 


KKF2.I -0.95 




615 


0.60:0.03 


F1>F2.I 0.48 




Inlvr-population 




















p 


42 


1 X: 11 15 


P>FI.t»l.35 




42 


1.00x0.13 


P 1 1.1 I 01 




ll\() 


Fl 


194 


11.85:0.05 


P -F2.1 4.50 




190 


0.8710.05 


P -F2.I 4,80 


mm* 




F2 


4.3.3 


11.49:0.0! 


II 1 2.1 6 0 5 




4m. 


04! 004 


11 1 2,t " 1)7 


*•* 




P 


46 


1) 85 Oil 


P 1 l.i o -0 




46 


0.69 0.09 


P 1 l.i -0 12 




IKixO 


1 1 


214 


II "" -1)115 


P 1 2.1 i S5 




20.3 


0 "1 0 05 


P 1 2.1 1 50 






F2 


446 


in,: oil! 


FI>F2.t=2.l9 




373 


0.54x0.03 


Fl>F2.t-2.86 






P 


44 


1.07x0.14 


P>F 1.1=1.91 




44 


1.02x0.10 


P>Fl.t=2.44 


* 


BGxB 


Fl 


207 


11 S 3 : 0.115 


P>F2.t=2.57 




215 


0.76x0.04 


P>F2.t=l.7l 






F2 


669 


0 "5 on: 


1 1 -2.1 1.91 




648 


0.82x0.03 


|| i:.: 0 9* 





O = Oak, B = Beech ,BG = Botanical garden. p<0.05 * p<0.0l **, p<0.00l *** 
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Table 9b. The FA I index differences between generations and type of crosses for sternopleural bristles in females (direct and reciprocal 
crosses); t test were used for testing differences in means of SB. 



FA1= IR-LI 






Direct 


Reciprocal 


cross 


generation 


\ 


nii 'iii ■ sl 
Mil .111-^1 




P 


N 


nit-aniS!-" 
nil .iii.^i 


t test 


P 


Inlra-pnpulalion 




















P 


16 




p<pj (^.0,i 3 




16 


0 S 6 - ( 1 Is 

V'.*7s'~ W. 1 o 


p.__ pi 1^,0.88 




i )\() 


HI 


6- 


(I (>4 II 06 


I 1 F2 t 0 21 




67 


(1 7 0-0 0(> 


P>F2 I - 0 86 






F2 


200 


0.59 -o n4 


Fl >F2.t=0.62 




167 


0.43^0.04 


Fl>F2.t=3.38 


*** 




I 1 


1 


0.85±O.3O 


P>FI.t=1.07 




1 ^ 


0.38-^-0.24 


P- : Fl.t=-0.97 




BxB 


Fl 


62 


0 63+0 07 


P>F2.t=2.03 


* 


55 


n ss- o ox 


P<F2.t=-0.62 






F2 


145 


0.45±0.05 


F 1 >F2.t- 1 .95 




138 


0.5 1 ^0.06 


Fl >F2. 1-0.62 






P 


49 


0.47-r0.10 


p<p] t=>1.99 




49 


0.53-rO.lO 


p<Fl t=-0.89 




H(,\K(, 


1 1 


208 


( 1 67 - ( 1 (U 


p<P2 (=. | 54 




21 3 


0 61 - ( 1 114 


p<P2 \=- 1 .2 1 






i : 


- 1 li 


II 64 -Oil' 


Fl>F2.t=0.53 




6N5 


0 66*0 03 


F1<F"> [=.n 86 

I 1 1 4» a l W . O x 1 




Intcr-population 




















p 


■42 


0.69=0.12 


P l i.i -1 93 




42 


0.74^0.09 


1' 1 1.1 -I.4S 




BxO 


Fl 


192 


().*>2 0.05 


i' i :.i i.76 




is: 


0.9U0.05 


P>F2.l-2.07 






F2 


516 


0.48^0.03 


F1>F2,K7.16 


*** 


493 


0.50^0.03 


F1>F2.I=6.59 


*** 




P 


46 


0.65^0.09 


P<Fl,t=-1.18 




46 


07S-0.I0 


P 1 1.1 11.0') 




BGxO 


Fl 


196 


u.^i 0.05 


P<F2,t--0.14 




185 


0 "^:0.O5 


P>F2,t=1.85 






F2 


527 


0.67±0.03 


I I 1-2.1 -I is 




482 


0.58=0.03 


F1>F2.I=3.19 


• * 




P 


44 


1.07^0.14 


P -Fl.t=l.91 




44 


0.84^0.11 


P=-Fl.t=0.07 




BGxB 


Fl 


207 


0.83^0.05 


P>F2.t-2.57 


* 


191 


0.83^0.05 


P>F2.I=0.20 






F2 


669 


0.-5-0.03 


Fl F2.i 1.32 




648 


O.S2:0.03 


PI F2.I 0.26 





O = Oak, B = Beech ,BG = Botanical garden. p<0.05 *, p<0.0l **, p<0.00l 



Table I Oa. The FA2 index differences between generations and type of crosses for sternopleural bristles in males (direct and reciprocal 
crosses); t test were used for testing differences in mean of SB. 



F\2=mcan R-l. mean , !( 1 : : 






Direct 


Reciprocal 


cross 


generation 


N 


nieaniSF 


ttest 


P 


N 


mean LSK 


(test 


P 


Inlra-populalion 




















P 


16 


(1 10111)4 


P I l.l 0.56 




16 


0.16:0,03 


PF2.I 2.84 




1 M I 


Fl 


7.1 


0.09:0.01 


P F2.1 1.09 




"1 


11.08 0.01 


P 1 2.1 3.66 


■ - - 




f: 


173 


0.07:0.01 


Fl -F2.1 1.05 




132 


0.1)7 -(1.01 


Fl 12.1 1.08 






p 


13 


0.09:0.01 


P Fl.t 1124 




1 5 


0.12*0.02 


P Fl.l 0.97 




BxB 


Fl 


61 


0.08*0.01 


p i ;.t o.5i 




58 


0.09:0.01 


P 12.1 I9| 






i : 


115 


0.07:0.01 


Fl F2.t 0.49 




1 Id 


0.06:0.01 


Fl 12.1=1.91 






p 


49 


0.11:0.02 


P>Fl.t=1.4() 




49 


0.07*0.02 


P- Fl.l=-I.32 




BGxBG 


n 


:i5 


0.09:0.01 


i' i :.i 1 4s 






0.101(1.01 


P 1 2.1 II 8') 






i : 




0.08:0.01 


i i i :.i -o ils 




615 


IMWi 0.01 


1 ! F2.1 1 1)6 




Intcr-population 




















p 


42 


0.18:0.03 


P>F1.I=I.79 




42 


(1 16 0.02 


P 1 l.l II 51 




lixl > 


Fl 


195 


0.14:0.01 


P>F2,t=5.41 




186 


(I 14:0.01 


P 1 2.1 4.67 


*** 




F2 


433 


D OS • 0.01 


Fl-F2.t-6.44 


*•• 


410 


0.07.(I0| 


1 1 1 2.1 7 5ii 


• ** 




P 


46 


O.I4±0.02 


P -Fl.t- 1.49 




46 


0. 1 1 ±0.02 


PI 1.1-0.01 




BGxO 


Fl 


215 


0.12*0.01 


PF2.t-2.54 




203 


0.11*0.01 


P>F2.t=l.88 






F2 


444 


0.10*0.01 


FI>F2.t 1.96 




372 


0.08*0.01 


F1>F2.1 3.17 


** 




P 


44 


0.12*0.02 


P- Fl.t -0.33 




44 


0.16*0.02 


P Fl.l 2 22 


* 


BGxB 


Fl 


208 


0.12*0.01 


P>F2.t 0.35 




215 


0.1 2 ±0.01 


P>F2.F2.60 


** 




F2 


582 


0.11*0.01 


Fl>F2,t-1.40 




582 


0.11 ±0.01 


Fl-F2.t0.74 





O = Oak population, B = Beech population, BG = Botanical garden population. p<0.05 *, p<0.0l**, p<0.00l*** 



Table 1 0b. The FA2 index differences between generations and type of crosses for sternopleural bristles in females (direct and reciprocal 
crosses); t test were used for testing differences in means of SB. 



FA2-mean R-I.l /mean |(R+I.V2| 






Direct 


Reciprocal 


cross 


"cneralion 


N 


mean±SE 


1 lest 


P 


N 


mean±SE 


ttest 


P 


Intra-population 




















P 


16 


0 111 11,111 


P<Fl.t=-O.I6 




16 


0.09*0.03 


P<Fl.t=-0.52 




OxO 


1 1 


62 


n in n in 


P>F2.t=0.47 




67 


ii 1 1 ii in 


P>F2.t=l.36 






F2 


199 


0.09*0.01 


Fl>F2.t=1.18 




167 


0.06*0.0 1 


Fl>F2.t=3.69 


*** 




P 


13 


(1.14-11.06 


P>Fl.t=1.25 




13 


0.07±0.()4 


P 1 l.l -0.61 




l!\B 


Fl 


1.2 


ll in mil 


P>F2.t-2.60 


• 


55 


0.09 0.01 


p 1 2.i -o ;o 






F2 


144 


0.06-0.01 


Fl>F2.t=2.41 


• 


138 


0.08*0.01 


Fl>F2.t=0.60 






P 


49 


O.0StO.02 


P<Fl.t=-l.43 




4'l 


0.1 1 ±0.02 


P>F 1.1=0.99 




Bl,\l!(, 


Fl 


207 


0 .10-0.01 


P<F2.t=-0.7l 




208 


0.10:0.01 


P>F2,t=l.21 






F2 


712 


0.09*0.01 


F1>F2.I"I,58 




686 


0.09*0.01 


F1>F2.I=0.37 




Inter-population 




















P 


42 


oil -11.02 


P<Fl.t=-l.94 




42 


0.12+0.01 


P<Fl.t=-1.54 




B\( l 


Fl 


192 


0.15x0.01 


P>F2.t=2.00 




180 


0.15±0.0I 


P 12. 1 1.91 






12 


518 


O.0S±O.0l 


Fl>F2.t=7.88 


... 


493 


0.08±0.0l 


F1>F2. 1=6.39 






!' 


46 


II 1 1 


P<Fl.t-l.27 




46 


0.13*0.02 


P -Fl.t -0.58 




BGxO 


1 1 


I9S 


0.1270.01 


PF2.t=0.25 




18" 


0.12i0.0l 


P 12. i 2 i' - 


• * 




F2 


528 


O.IO±O.OI 


F1>F2.I=3.08 




484 


0.09±0.0I 


F1>F2. 1=3.68 


*** 




P 


44 


0.I6±0.02 


p-Fl.t-l.OI 




44 


0.13±0.02 


P<Fl.t=-0.41 




BGxB 


1 1 


202 


0.I4±0.0I 


P=-F2.t=3.10 




190 


0.14*0.01 


P -F2. t 0.61 






F2 


673 


0.1U0.0I 


Fl>F2.t-3.58 




648 


0.12*0.01 


FI>F2. 1=1.97 





O = Oak population, B = Beech population, BG = Botanical garden population. p<0.05 *, p<0.0l**, p<0.00M 
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